On the dynamics of spin systems in the Landau- Lifshitz theory 



In the framework of the Landau-Lifshitz equations without any dissipation (an approxima- 
tion which may also be helpful for finite but weak Gilbert damping) , with all interactions 
included (specifically always with exchange and magnetic dipole interactions), for general 
ground states, geometries and domain structures, and for many types of effective fields (in- 
cluding contributions from spin currents etc.) the dynamics of the spin-precession around 
this ground state is considered. 

At first the precession is treated in the linear approximation. For the (small amplitude) 
eigenmodes of the precession one has a 'rule of geometric mean' for the eigenfrequencies. 
For the eigenmodes pseudo-orthogonality relations are obtained, which reflect the gy- 
rotropic and elliptic character of the spin precession and differ from those known from 
the Schrodinger equation. Moreover, in contrast to quantum mechanics, they are defined 
locally, not through a global average. Thus pseudo-orthogonality statements are valid ev- 
erywhere (for example simultaneously both in the outer region and also in the core region 
of a magnetic vortex line) . 

Then also some aspects of the nonlinear mode coupling with emphasis on 'confluence' 
and 'splitting' processes of elementary magnetic spin-wave excitations are considered. At 
the same time these processes contribute to the Gilbert damping. 

There are thus essential differences to quantum mechanics, although at a first glance 
one discovers many similarities. From the results one may also get insights of why these 
systems are so complex that (although the essential quantities depend only on the lo- 
cal values of the partially long-ranged effective magnetic fields) practically only detailed 
experiments and computer simulations make sense. 

PACS numbers : 75. (Magnetic properties); 75.40 GB (Dynamic properties) 
Keywords: Landau-Lifshitz Equs., Dynamics of Magnetic Systems 



1 Introduction 

Recently, for possible applications in spintronics, the dynamic behaviour of ultra- 
thin planar magnetic nano-structures is strongly studied ([1,2,..., 9], and and references 
therein). For example, one considers circular structures with nontrivial ground states 
(notably with vortex states). 

Theoretically, on the spatial scale of typically ten (or more) nanometers (nm) and 
on the time scale of ten (or more) picoseconds (ps) the Landau-Lifshitz equations, |10j . 
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with small phenomeno logical 'Gilbert damping', is very reliable and has found 
many applications in present-day simulation software (for example, in the well-known 
OOMMF, [12], or LLG programs, [13J ) . It is 'general folklore' that there are consid- 
erable similarities between the Landau-Lifshitz-equation without Gilbert damping and 
the quantum mechanical Schrodinger equation. This would be advantageous, since for 
weak Gilbert damping one could then exploit a lot of quantum-mechanical relations. 

Unfortunately, however, difficulties arise, since (as we will see) there are considerable 
differences from the Schrodinger case. 

The present communication intends to make these points clearer; it also serves more 
general purposes, for example in the intention to go beyond numerical simulation as 
far as possible. 



2 Basic theory 

In the framework of the Landau-Lifshitz theory, [10], the local value of the magnetic 
polarisation J(f , t) (as a function of position r and time t) is described by the following 



where J s is the saturation magnitude and a(f, t) the direction of the magnetic polari- 



This fact is important throughout the paper. Essentially due to this constraint 
on a(f,t) the Landau-Lifshitz theory has the same topological exitations (Bloch lines, 
Bloch point singularities, vortex excitations, and all kinds of 'topological constrictions') 
as in field-theory the so-called 'non-linear sigma- model', [H], has; or as in elementary 
particle physics the corresponding phenomenological 'Nambu and Jona-Lasinio model' 
has, [15] . However, since the magnetic dipole field (see below) is far-ranged, the Landau- 
Lifshitz theory is in principle even more complicated than the above-mentioned field 
theories (which usually do not treat interactions corresponding to the dipole-dipole 
case) . 

In any case, the role of the topological constraints in applications (for example, in 
the wall-motion of magnetic bubble- domains, [16] , or in the change of the polarization 
of the core of a magnetic vortex-line, [TF]) should not be underestimated. 

However the present paper does not concentrate on topological aspects ! 
Rather, in the linear approximation (neglecting the damping unless otherwise stated) 
we derive a 'rule of geometric mean' for the eigenfrequencies and a pseudo-orthogonality 
relation for the eigenmodes of the precession around the nontrivial ground state. The 
relations correspond to a gyrotropic and elliptic motion of the magnetization; they differ 
from those known from quantum mechanics. In particular they are defined locally. 

Nonlinear terms are treated in the present paper rather phenomenologically. We 
only give some schematic arguments for damping processes based on 'confluence' and 
'splitting' in connection with the so-called mode-mode coupling of excitations in vortex 
structures. The results are in favour of the experiments of [2] and [3J. 



ansatz: 





sation vector J. The vector a is thus of unit length. 
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3 Landau-Lifshitz equations 



In the absence of any dissipation (see above), on the above-mentioned spatial and 
temporal scales, the equation of motion for the vector field a{f,t) is the undamped 
Landau-Lifshitz equation, [10], which can be written as 



= -j . a(r, t ) x H eS (r,t) , (2) 

where the positive quantity 7 = g • is the so-called Lande factor g (« 2) multiplied 
by the gyromagnetic ratio; e (= — |e|) is the electron's charge (a negative quantity), m e 
its mass, and /xo the permeability of the vacuum; H e s(r,t) is an effective field, namely 
the sum 

H cS := + H di ^ + £ ^ (J^j ■ (J^j - grad s F (a) + i? thcrm + W . (3) 

In ([3]), H cxt is the external magnetic field; the are the so-called exchange con- 
stants of the crystal; F (a) is the anisotropy energy, representing the easy axes of 
the magnetization; jj therm is a fluctuating thermal field oc VT, where T is the Kelvin 
temperature; this field is neglected together with the friction, to which it is intrinsi- 
cally related. Finally, Hi is an effective field proportional to the density of the electric 
current, j, or to other sources of torques acting on the magnetization. For stationary 
currents all these contributions to H e s do not depend explicitly on the time nor on the 
time derivative, [19]. Often they all can be neglected with respect to the second and/or 
third term on the r.h.s. of equation (J3]). 

The remaining second term on the r.h.s. of ([3]), H dipole , is the cumbersome quantity 
in micromagnetism, namely the magnetic dipole field. In contrast, the third term 
represents the exchange interactions, which at short wavelengths are much stronger, 
such that practically they alone determine the critical temperature T c of the ferromagnet 
considered. But for applications at room temperatures and below, at distances which 
are much larger than the atomic lattice constants (typically 0.2 nm) the dipole fields 
are most important. 

Therefore even in the /im and nm ranges, as mentioned, the exchange interactions 
can often be neglected (we don't neglect them here!), whereas the field H dipole , which 
is, for example, mainly responsible for the domain structure of a ferromagnet, can be 
derived from the relation H dipole = — grad$ m . 

The magnetostatic potential $ m is determined as follows: 

(i) either from the magnetic moment representation 

(here the integration volume V is that of the magnetic sample), 
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(ii) or from the equivalent representation of the fictitious magnetic charges : 



® fffv dV - div , J i f '-i> + #„ m f'y;' . ( 5) 

' ' 47r/io |r — r | i/ 47TyUo |r — r | 



Here dV is the (oriented) boundary of V; the vector n is the its outer normal ; 
the surface measure is d 2 A. (Retardation effects do not play a role, unless the 
frequencies / (= ^-) are larger than ^- (= 3 • 10 15 Hz for Ax = 100 nm ).) 

In any case, the local value H dipole (f,t) of the dipole field depends on the magnetiza- 
tion distribution of the whole sample; i.e., in the context of the Landau- Lifshitz equation 
the dipole field is a long-ranged, but local quantity, which sounds as a contradiction in 
itself, but becomes nonetheless essential below. 

Unfortunately, the magnitude of this field can also be very large, namely at surfaces 
typically as large as j^, which corresponds to ~ 2.2 Tesla in Fe, ~ 1.7 Tesla in Co and 
~ 1.0 Tesla in permalloy So for applications the magnetic dipole field is as important 
as the exchange. 



4 Basic relations for the dynamics 
4.1 General equations 

In the following we write a = a + j3, where a (r) is the static ground state, correspond- 
ing, e.g., to a non-trivial domain structure of the sample, whereas the vector /3(r,t) 
describes the dynamics of the sample. Here we use the linear approximation with re- 
spect to the ground state (i.e., we assume that \/3\ is <C 1, while oioif) ■ /3(f,t) = and 
c?o(r) = 1). As a consequence in linear approximation also a 2 (r,t) = 1, as desired. 

Furthermore we don't exclude topological singularities, e.g., Bloch points or vortex 
lines, but assume that these are exclusively contained in the statics of the system (i.e., 
in a (r)), whereas j3(r,t) is topologically trivial (but see [20J). 

In the linear approximation one gets the following equation of motion : 

- 42 = 5 (f) x Hi + P(r, t) x H2 . (6) 

Here we have used the reduced time t := t • 7 (in the following we replace 7 by 1, unless 
otherwise stated). 

But H^ff {f) has the same direction as (r) (the reader is reminded to the relation 
H c q = 0). Therefore the final term in (Q can be replaced by a term with a 
Lagrange parameter function A a ° (r) : 

- ^ = 5 (f) x (i?i(r,t) - A^(f)/3(f,t)) . (7) 

This is a linear relation in f3. Generally only the local and instantaneous values of the 

effective field H^ s (and of 0) are involved. It should also be noted that d7j) does not 
explicitly depend on t, see (@J or <^.) 
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Now we assume the existence of local cartesian coordinates such that the local 3- 
axis is ao(f), whereas e*i(r) and e*2(r) are perpendicular to c?o(r) and to each other, but 
otherwise not yet fixed : Then in linear approximation the following equation for the 
vector ftf, t) (which also must be perpendicular to c?o) results : 

-^ t tPx(?,t)e x {f) + p 2 {r,t)e 2 (?)) =K(f) (ft (ft) e x {f) + ft(f, t) e 2 (f)) . (8) 

Here K (r) is an antisymmetric real 2 x 2- matrix, i.e., with imaginary eigenvalues 
(see below), such that (after suitable rotation of the axes perpendicular to a (r)) the 
equations for ft(f, t) and ft(r, t) become of the simple form 

-^ = +h e 2 ff (?)p2 (9) 

and 

Pi, (10) 



dt 

where hi* (r) and h e 2 ^ (r) are extremal off-diagonal values of the 2 x 2-matrix K if) , 
see below, which we assume to be non-negative, unless otherwise stated. 

The equations (jSJ) and ( fTUI) describe an elliptic spin precession around the local 
equilibrium spin direction do(r), and the quantities h\ and h 2 are local effective 
fields corresponding to the lengths of the principal axes of the local precession ellipse. 

In the following we skip the superscripts of hi and h 2 . 



4.2 The 'rule of the geometric mean' for the eigenfrequencies 

We now make the ansatz ft(f, t) = ft (r) cos(u;t) and ft(f, t) = ±f3%±(f) sin(a;t), 
where all quantities are real, to is positive, unless otherwise stated, and both polari- 
sations enter (the upper (rsp. lower) sign in the ansatz for ft correspond to mathe- 
matically positive (+) (rsp. negative, (-)) circularly polarised precession; typically one 
has a superposition of both polarisations, which results in the above-mentioned elliptic 
precession. If the external field iff xt (f) and/or the exchange field dominate the local 
effective field, see below, then one has purely positive (for positive i^3 Xt ) (rsp. negative, 
for negative -£^ xt ) circular polarisation.). 

In any case, by the product of (Q and (fTOj) one gets finally the eigenfrequency 
equation 

M 2 ■ sin w v t ■ (± cos u v t) ■ (3{ 0) (f) • ^ 0) (f) 

= 7 2 • h^\f,t) ■ h ( f\r,t) ■ cos u„t ■ (isinc^t) • ^ 0) (f) • p{°\f) , (11) 

where the gyromagnetic ratio 7 has been restored. For the eigenfrequencies this 
yields the (well-known, see [2TJ) rule of the geometric mean (the index v counts the 
eigenstates) : 

u u = ^ h^\r) ■ h^\r) . (12) 

Here one should note that the l.h.s. of the result does not depend on f, in contrast 
to the r.h.s. In fact in linear approximation, in the Landau-Lifshitz theory, one can 
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perform any position-averaging of the time-dependent signal (for example, over all spins, 
but usually with different weights at different positions) : the eigenfrequency spectrum 
will always be the same. (This is useful for applications, [3], and it also allows that 
even from a distant position there may be contributions to H{? s (f,t).) 

The relation ffl2|) is actually at the same time more complex and also more simple 
than in quantum mechanics: on the one hand one has two functions, hi(r) and h,2(r), 
instead of only one in quantum mechanics (i.e., the potential energy V(r)), whereas in 
quantum mechanics, on the other hand, not all space-averages, but only the full-space 
average with the eigenfunction, / d 3 r ip*(r) Ti{r)%l) v {r), corresponds to Tluj u (in the last 
two equations all quantities have their usual meaning). 

If we multiply the canonical representation of the 2 x 2-matrix K (r), 

s«=(-4, : h f >. («) 

(where the 1- and 2-directions are perpendicular to a^°'(r) and extremal (see below)) 
by the imaginary number j, then equation (JH]) becomes similar to a Schrodinger equa- 
tion with two-component vectors, (e.g., with spin-orbit interaction), although in our 
case the quantities ^4-^ (for j = 1,2) are independent imaginary functions, not nec- 
essarily fulfilling the identity (^M)* = _^il£2 (j_ e _ ; hermiticity) : only for h 2 = hi 



[circular precession) one would obtain K . ^ = a y ■ li2(r) (i.e., a 2 x 2-Schrodinger equa- 

, -i 
+i , 



tion involving the (hermitian) Pauli matrix a y — ( ( . ' n ) ). Instead, for elliptic 



precession one obtains K ^ r - > = b y ■ h2( -^>+ hl ^ _)_ fr x . h i(r) ^ 2 (r) _ • ^ ^ _ f > J ^ 

Thus in the elliptic case, the matrix ^p- is no longer hermitian, but still pseudo- 
hermitian^, in the sense of the just-mentioned identity; see [22]. The hermiticity would 
be lost from the beginning, if the Gilbert damping were included 



4.3 Pseudo-orthogonality of the eigenmodes 

Now we complexify the eigenmodes (3 U as follows : Starting from (3^ (r , t) = (r) • 
cos(ccvt) + $\r) ■ sin(^t) we write cos(^t) = Tie (e~ iuJl/t ) and sm^tu^t) = Tie (e~ iWl/ 'i). 
In this way (by omitting the real part Tie) one obtains the primed quantity (3^ ' (r, t) : = 
Pi \r) ■ e~ lWvt + (3% (f) ■ e~ lWvt \ (i.e., oc e~ lWvt ). By expansion with the complete set 
of eigenmodes this can be extended to general vectors /3(r,t). One also sees by direct 
inspection that an eigenvector (3' cx e~ lut of the equation —d(3'/dt =K (3' with the 
eigenfrequency uo u is obtained from the compatibility relation f3 2 ' = ih^^h 1 ^ 2 (3[. 

In the Landau-Lifshitz formalism (denoted by the symbol LL) the pseudo-scalar 
product of two vectors (3 can now be defined as follows (i.e., locally, and with the 

J It can be shown that this is equivalent to the statement that the diagonalisation can only be performed 
with a non-trivial Bogoliubov-Valatin transformation, see |21j . 
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primed quantities on the r.h.s.) : 

{fr\rt)\fr\r,t)) LL := (3{ vy (r,ty ■ h 2 (r)& r (r,t) - $ y \f,t)* ■ h t (f) \r,t) , 

(14) 

where the *-symbol denotes the hermitian-conjugate quantity. 
In this way one gets (3* oc e +1£iV *, and with (TX5T) 

{fr\rt)\fr\ft)) LL := fr y (f ,t)* ■ (K (r)fr y (f,t)) . (15) 

This generalizes a relation of W. Fuller Brown, Jr., [TS], which was restricted to homo- 
geneous ground-states. 

The main result of the present formalism is of course the following statement : 

Two eigenmodes fr^ rsp. fr> with different eigenfrequencies, ^ w^, are ev- 
erywhere pseudo-orthogonal. 

E.g., the so-called radial and azimuthal eigenmodes, respectively, of pQ are every- 
where pseudo-orthogonal, both in the outer region and also in the core region of a 
magnetic vortex line; although in the outer region a = ±e v , Pi = ±/3 z , (3 2 = (3 r , 
whereas in the core region a w ±e 2 , (3\ ~ ±/3 r , /3 2 = Z^. 

This pseudo-orthogonality, which reflects the gyrotropic and elliptic motion of the 
spins, is also valid for any spatial average, see |23j. 

Similarly to the formalism of quantum mechanics, but with the additional freedom 
of locality, one should thus be able to estimate the eigenfrequencies from variational 
calculations, for example by exploiting equation ([12]) for points with suitable properties. 

The proof of the above-mentioned main result is straightforward (here the super- 
script t means transposed, and it should be noted that the primed quantities are used) : 

P {u) '(r,t)* ■ (K (f)/^'(f,t)) = -iu^y* 

= (K' (f)fr yt )-fr y = (-K(f)fr yt )-fr y = -iuJ {uyt -fr y , (16) 

i.e., one minus sign is produced by the transposition, and a second one by the hermitian- 
conjugation. As a consequence, the states are pseudo-orthogonal if the eigenfrequencies 
differ. 

Another essential point of the result is not the above-mentioned strange 'preces- 
sional' form of the pseudo-scalar product, but rather the fact that in contrast to quan- 
tum mechanics, which is nonlocal (see also [21]), here the local values of classical 
effective quantities count. 

We can again apply the pseudo-orthogonality relations to vortex states. Here, suffi- 
ciently far from the vortex core, one has a(r,t) already mentioned. For h\(r) 
and h2{r) one should see the remarks below. 

5 Further remarks on the relation to quantum mechanics 

Some statements on the relation to quantum mechanics have already been presented. 
The following is more subtle and at the same time more basic : the Landau-Lifshitz 
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equations form a nonlinear set of classical integro-differential equations. The final quan- 
tities f3i(r,t) and ^(r, t) are real (only if one wants to introduce the above-mentioned 
pseudo-orthogonality one is apparently forced to make the ansatz oc e~ luJt ) and in linear 
approximation they can be taken from a real Hilbert space. In contrast, in quantum me- 
chanics, which is totally linear, these observables are represented by hermitian operators 
acting on the elements of a complex Hilbert space. Although the Landau-Lifshitz equa- 
tions can be derived from a quantum mechanical model (e.g., the Heisenberg model) 
this involves some kind of molecular-field approximation (e.g., products of expectation 
values, (A)(B), instead of expectation values of products of operators, (AB)). One 
should of course also consider metallic ferromagnets, which in quantum mechanics are 
better described by an itinerant model instead of the simpler Heisenberg one. But in 
any case both correspond to the same phenomenological Landau-Lifshitz theory. In 
the 'itinerant case' the atomistic derivation of the phenomenological theory is more 
delicate (for the exchange field the problems have been overcome by Korenman and 
coworkers, [25]). Therefore, although it was shown above that the equation of mo- 
tion for the classical vector f3(f, t) in linear approximation looks very similar to the 
Schodinger equation, the correspondence to the Bogoliubov-Valatin approach (see [2"Tj ) 
is actually only true if one performs a quantum mechanical approach from which the 
Landau-Lifshitz equations can be derived. 

A quantum mechanical model is also the basis of the following section. 

6 Nonlinearity : mode-coupling, confluence and splitting 

Now the third-order processes, which come into play for higher oscillation amplitudes, 
are treated rather schematically in an approximation corresponding in quantum me- 
chanics to the Bogoliubov-Valatin formalism : if bf and b\ are the well-known creation 
and destruction operators for (Bogoliubov-Valatin) quasi-particle mode i, then in the 
next order the Hamilton operator contains a sum (or integral) of additional terms of 
the form 8H := h- • 63(61" b^) + h + ■ (b 2 bi)bf , while h- and it's hermitian-conjugate 
h + are complex amplitudes of physical dimension 'energy'. (In the following it is not 
necessary to calculate these quantities from a(f, t), although just this is implicitly done 
in the computer simulations.) 

According to Fermi's 'Golden Rule' the relaxation times for the corresponding processes 
derive from integrals involving the expression 

(-L) := ^ . l/fcpf . 5(E 3 -E 1 -E 2 ), (17) 

T T II 

where splitting and confluence processes correspond to the — or + cases, respectively. 
(In principle this is also not new, see [29].) Of course one has \h + \ 2 = |/i+| 2 = \h-\ 2 . 
However, the populations are different (see below). 

In the above h T factor and the (E3 — ...) argument, respectively, one has both a 
spatial and also a temporal constraint corresponding to the mode form and frequency, 
symbolized by '3 — > 1 + 2' ('splitting') and '1 + 2 — > 3' ('confluence'), with the en- 
ergy corresponding to the frequency (through the usual relation AE = HAuj). These 
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contribute to the 'Stokes' and 'Antistokes' populations, respectively. At low tempera- 
tures the Stokes population ('splitting') is frozen, since excitations to be split do not 
exist. Therefore through the temperature dependence one can check the essentials of 
the mode-coupling theory. 

In fact, in this way one can simultaneously understand parts of the above-mentioned 
phenomenological generalisation of the Landau-Lifshitz equation by Gilbert's damping 
term, [11] . Moreover, (only) if after integration the typical relaxation processes do not 
distinguish a particular axis, the Gilbert damping is isotropic per ansatz. 

7 Application to vortex states (radial and azimuthal modes) 

In the following we consider in more detail a flat circular disk with radius R of the order 
of 1 fx and thickness th <C R. Outside a central region of radius O(Xe) '■= ®(J~^jz^) 
(\e = ^/ jP ~ 5.7 A for Permalloy) the ground state a of the disc is the vortex 
state (i.e., do(r) = e v , where e v is the azimuthal unit vector). The precession around 
the vortex state is strongly elliptic, almost in-plane. The transverse fields entering 
the rule of the geometric mean are partially known : The out-of-plane component is 
h\ = f 1 + (if no external field is applied). In magnitude, this is a very large 

field (/io|/ii| ~ 1 T for permalloy) driving the spins backwards into in-plane direction, 
and practically independent of the mode index v (but see [26]). This is already an 
exceptional case: to keep /3(r, t) stationary (i.e., u real), also h 2 must be negative. 

Moreover, in contrast to the out-of plane value h\, the in-plane component hi de- 
pends strongly on v. If the spatial pattern of the mode corresponds to a plane wave 
with (in-plane) polarisation parallel to the magnetization, i.e., k ~ cto(r) (or for an 
azimuthal excitation of a vortex configuration, ~ e v ) then h 2 almost vanishes, since 
there are no (effective) magnetic charges generated in the bulk (again: but see [26J 
concerning the so-called backward modes); whereas in case of an (in-plane) variation 
corresponding to a propagation k _L 3o(r) (or for a radial excitation) the effective 
magnetic field h2{r) is maximal. 

As already mentioned, the experiments of [3] show the existence of radial and az- 
imuthal modes, respectively. In agreement with the above statements the eigenfrequen- 
cies of the azimuthal modes are relatively low. In fact, for the lowest modes they are 
less than half as high as that of the mode corresponding to the main radial peak. Thus, 
energetically, there can be confluence processes (at low temperatures and at room tem- 
perature) and splitting processes (only at room temperature, where they are equally 
probable as the confluence processes), involving for example a splitting of a radial main 
mode into two azimuthal modes. (Of course this tentative suggestion of the quenching 
of these processes at low temperatures has to be carefully checked experimentally). 

8 Conclusion 

In this communication some results (mostly not new, but nonetheless probably useful) 
for the dynamical properties contained in the Landau-Lifshitz equation for ferromag- 
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netic systems with general ground state have been demonstrated. It was shown that 
the linear approximation of these equations has strong similarities but also strong dif- 
ferences compared with the quantum mechanical Schrodinger equation. The essential 
difference is 

1) perhaps not the above rule of the geometric mean, 

2) also not the complications arising from the magnetic dipole interaction, 

3) but the simple fact that quantum mechanics is nonlocal (e.g., the state function 

comprises all values ip(r) throughout the space), whereas the Landau- Lifshitz 
theory is a classical local theory (i.e., the local values H c s(r) of the effective fields 
count, although the sources of these fields may be at a very large distance). 

Generally one has got relations which are more complicated than in the quantum me- 
chanical case. We mention (i) again the 'rule of the geometric mean' for the eigen- 
frequencies and (ii) a local pseudo-orthogonality of the eigenmodes, reflecting the gy- 
rotropic and elliptic properties of the precession. Moreover, the Gilbert damping was 
partially traced back to microscopic three-( or more-)particle processes which corre- 
spond to confluence or splitting of elementary magnetic excitations as seen in recent 
experiments, [3]. 

Essential results depend on the local properties of the magnetization structure (for 
example, a domain structure). This can be at the same time both an advantage and also 
a disadvantage of the formalism; i.e., in addition to the known problems of calculating 
the magnetic dipole fields one sees at this place explicitly why these systems are so 
complex that one can hardly avoid a comparison with detailed numerical simulations. 
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